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Electron spin relaxation in graphene on a substrate is investigated from the microscopic kinetic 
spin Bloch equation approach. All the relevant scatterings, such as the electron-impurity, electron- 
acoustic-phonon, electron-optical-phonon, electron-remote-interfacial-phonon, as well as electron- 
electron Coulomb scatterings, are explicitly included. Our study concentrates on clean intrinsic 
graphene, where the spin-orbit coupling from the adatoms can be neglected. We discuss the effect of 
the electron-electron Coulomb interaction on spin relaxation under various conditions. It is shown 
that the electron-electron Coulomb scattering plays an important role in spin relaxation at high 
temperature. We also find a significant increase of the spin relaxation time for high spin polarization 
even at room temperature, which is due to the Coulomb Hartree-Fock contribution-induced efi'ective 
longitudinal magnetic field. It is also discovered that the spin relaxation time increases with the in- 
plane electric field due to the hot-electron effect, which is different from the non-monotonic behavior 
in semiconductors. Moreover, we show that the electron-electron Coulomb scattering in graphene is 
not strong enough to establish the steady-state hot-electron distribution widely used in the literature 
and an alternative approximate one is proposed based on our computation. 

PACS numbers: 72.25.Rb, 71.10.-w, 73.61.Wp 



I. INTRODUCTION 

Graphene, as a strictly two-dimensional material, has 
revealed a cornucopia of new physics and potential appli- 
cations, and thus has attracted much attention in recent 
yearsii This material is also important for spintronics 
since the spin relaxation time (SRT) in intrinsic graphene 
is expected to be very long. The underlying reason is the 
low hyperfine interaction of the spin with the carbon nu- 
clei (natural carbon only contains 1% ^^C with spin) and 
the weak spin-orbit coupling (SOC) due to the low atomic 
number.—""— 

The study of spin relaxation in graphene is still in the 
initial stage. Some investigations have been performed 
on the spin relaxation due to the D'yakonov-Perel' (DP) 
mechanisrr^ in graphene i^""— This relaxation mechanism 
is from the joint effects of the momentum scattering 
and the momentum-dependent spin-orbit field (inhomo- 
geneous broadenin g^^d'^ ). However, the previous inves- 
tigations are all in the framework of single-particle ap- 
proach, thus the electron-electron Coulomb scattering, 
which has been demonstrated to be very important for 
the spin relaxation in bulk and low-dimensional semicon- 
ductor systems,—""— was not included. Understanding 
the effect of the electron-electron Coulomb scattering on 
spin relaxation in graphene is an important problem. In 
addition, the Coulomb Hartree-Fock (HF) term acts as 
an effective longitudinal magnetic field, and hence can 
increase the SRT by more than one order of magnitude 
for high initial spin polarization in semiconductors at low 
temperature:— ""—i^S Whether it is still valid in graphene 
remains unchecked. Also, in semiconductors the spin re- 
laxation can be effectively manipulated by the high in- 



plane electric field. -'^^^^^ How the in-plane electric field 
affects the spin relaxation in graphene is also unclear. In 
the present paper, we investigate the spin relaxation in 
graphene from the microscopic kinetic spin Bloch equa- 
tion (KSBE) approach, ^'^ which has achieved much suc- 
cess in the study of the spin dynamics in semiconduc- 
tors. Via this approach, we can explicitly include all 
the relevant scatterings, especially the electron-electron 
Coulomb scattering, and understand the physics raised 
above. 

It is also noted that there is a significant discrepancy 
between the existing theories and the recent spin trans- 
port experimentsi^""— These experiments reported the 
SRTs of only about 150 ps, at least one order of mag- 
nitude shorter than the lowest value obtained in the 
theory.—""— This suggests that the SRTs obtained in the 
recent experiments are likely to be limited by an extrin- 
sic mechanisms, e.g., the local spin-orbit field from the 
adatomsjiSiiii^ In this paper we concentrate on the rela- 
tively cleaner graphene samples by choosing low impurity 
densities which give a mobility higher than the values in 
the latest experiments, so that the effect of the adatoms 
can be neglected. 

This paper is organized as follows. In Sec. II, we 
present the model and introduce the KSBEs. Then in 
Sec. IIIA, we discuss the effect of the electron-electron 
Coulomb interaction on spin relaxation at various tem- 
perature, electron density, initial spin polarization and 
in-plane electric field. We discuss the hot-electron distri- 
bution function in the steady state in Sec. IIIB. Finally, 
we summarize in Sec. IV. 
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II. MODEL AND KSBES 

We start our investigation from a graphene layer on a 
Si02 substrate. The z-axis is set perpendicular to the 
graphene plane. A uniform electric field E|| and a uni- 
form magnetic field B are applied along the x- and y-axes 
respectively (the Voigt configuration). Without the SOC 
and the external field, the band structure of graphene 
near the K and K' points can be described by the effec- 
tive Hamiltonian (ft = 1)^ 



(Tyky) 



(1) 



Here fi — 1(— 1) for K{K') valleys; wf is the Fermi veloc- 
ity; k represents the two-dimensional wave vector relative 
to K(K') points; cr is the Pauli matrix in the pseudospin 
space formed by the A and B sublattices of the honey- 
comb lattice. The eigenvalues of Hq are e^i/k = i^WF|k| 
with v = 1(— 1) for electron (hole) band. The corre- 
sponding eigenstates are ■i/'k" ~ Xj ^fl^iive""^^^^ , V)^ with 
0k representing the polar angle of k. We introduce an or- 
thogonal and complete basis set = tp^'^ (8) Xs with Xs 
being the eigenstate of the spin Pauli matrix Sz- In this 
basis set, the total Hamiltonian including the SOC can 
be written asii*^ 

HeS = ^ [(eAii^k - Ai - eE|| • IVjSss' 

+ (.9MbB + lyn^) ■ Sss'K'JcZ + Hint- (2) 



where R — (x, y) is the electron position, c- 



') is 



the annihilation (creation) operator of the state ^'j^^ . e is 
the electron charge (e > 0). g = 2 is the effective Lande 
factor. The intrinsic SOC coefficient Ai = 0.012 meV is 
from the recent first-principle calculation.'' From Eq. 
it is seen that the intrinsic SOC only induces a shift of the 
energy spectrum of graphene. fly^ denotes the effective 
magnetic field due to the Rashba SOC, which reads^i 



J^k = cir{— sin (?k, cos(?k, 0) 



(3) 



with ap, = C^z- The recent first-principle calculation 
gives C = 0.005 meV-nm/VJ The longitudinal electric 
field Ez originates from the gate voltage and chemical 
doping. Here we choose a typical value in experiment 



Ez = 300 kV/cm^iH It is noted that fly, depends on the 
direction of k only, but is independent on the magnitude 
of k. Therefore the inhomogeneous broadening induced 
by the SOC does not change with the variation of temper- 
ature and electron density. This makes the temperature 
and electron-density dependences of the SRT in graphene 
very different from those in semiconductors The inter- 
action Hamiltonian Hi^t consists of the electron-impurity, 
electron-phonon as well as electron-electron Coulomb in- 
teractions. Their expressions are given in the appendix. 
In the derivation of Eq. |e^i/k| ^ or + Ai is as- 
sumed and thus the terms between the electron and hole 
bands are neglected. This approximation is valid when 
the Fermi energy Ep is much larger than 0.03 meVr-i^ 
which is usually fulfilled in gated or doped graphene. In 
our calculation, we restrict ourselves to the n-doped case 
(i.e., Ef > kBT). 

By using the nonequilibrium Green function method)^ 
the KSBEs can be constructed as:- - 



coh 



9t/5^k 



drift 



(4) 



where p^k represent the density matrices of electron with 
the relative momentum k in valley /i, whose diagonal 



terms p^k,. 



fnks (s = ±^) represent the electron 



distribution functions and off-diagonal ones p^i^i_i ~ 
p*. _i 1 describe the spin coherence. The coherent term 
is given by 



9t/5;,klcoh = + rJk) • s 



VHF 
^A'k ' 



(5) 



in which [A, B] = AB - BA is the commutator; = 
~ Sk' ^k'^kk'P/jk' is the effective magnetic field from 
the Coulomb HF contribution.^^ The drift term can be 
written asi^ 



dtPt^ik 



drift 



= eEii • V 



kP/jk, 



(6) 



The scattering term 9t/5/ik|scat consists of the electron- 
impurity, electron-acoustic (AC)-phonon, electron- 
optical (OP)-phonon, electron-remote-interfacial (RI)- 
phonon as well as electron-electron Coulomb scatterings. 
These scattering terms are 



dtp.kl, = -nN, J2 Kk'P^kk'<5(epk' - e^k) (/5>k'P<k - P<k'/5^k) + H.c 



(7) 



dtPf,k\ 



E \M^k,,'k'\'Sis,^k' £pk ± ^^kO [Ntk-k'P^'P^k - ^A^k-k'/5<'k'P>k) + H.c, (8) 

/x'k'A,± 

— TT E l^k!k' P^kk'^k"k"-k+k"5(£A'k'— ^/ik+^At'k"— £/i'k"-k+k') 
fi'k'k" 



Tr ( P'i:'k"-k+k'Pl'k" ) P^k'/'^k 



Tr (/5>.k"-k-(-k'PM'k") P^k'P^k 



H.c. 



(9) 
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In these equations, 
1 - Pnu. P. 7,, = p^k; 

iVAq 



= lk 



P^k, P^k 

trum; N^^ - 



denotes the phonon energy spec- 
it 5 with Nxq representing the 



phonon number at lattice temperature. The form factor 
<ial' = ^[1 + cos(0k - ^k-q)]. The ma- 

A 

^k,^'k' 



/kk-q = IV'k^Vkiql - 2 



trix elements C/^k'' ^k' ^^'i -^^^k.^'k' ^'"^ given in the 



appendix. 



III. RESULTS 

The KSBEs with all the scatterings explicitly in- 
cluded can be solved self-consistently following the nu- 
merical scheme similar to that in semiconductors, de- 
tailed in the appendix of Ref. IH. Then one can ob- 
tain the temporal evolution of the electron density ma- 
trix. The SRT Ti and ensemble spin dephasing time T2 
can be determined from the slopes of the envelopes of 

AiV(i) = E;.k(/^ki -/;.k-i) and Pit) 
respectively] ^^1'^^ Since the Rashba spin-orbit field is in 
the graphene plane, these two SRTs satisfy Tj* — 2Ti 
in the case without magnetic field. In the presence of 
an in-plane magnetic field, the spin relaxation becomes 
isotropic, i.e., Ti = T* = 2/{Ti{B = 0)-^ +T*{B = 
0)-i) = |ri(B = 0)M. In the following, we only show 
the SRT T = Ti. Unless otherwise specified, we choose 
initial spin polarization P = 1 %, electron density N^. = 
7 X 10^^ cm~^ {Ep = 100 meV), external magnetic field 
B = and electric field E|| = 0. The effective impurity 
density is chosen to be A'^ = 2 x 10^^ cm~^. The corre- 
sponding mobility is /i = 3 x 10** cm^/V-s at 100 K, which 
is of the same order of magnitude as those reported in 
the experiments but one order of magnitude higher than 
those in the recent spin transport experiments;^^— 
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FIG. 1: (Color online) Total SRT (a) and the contributions 
from the electron-impurity (•), electron-electron Coulomb 
(□), electron-AC-phonon (o), electron-RI-phonon (t) as well 
as electron-OP-phonon ( x ) scatterings as function of temper- 
ature T. 
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FIG. 2: (Color online) Total SRT (a) and the contributions 
from the electron-impurity (•), electron-electron Coulomb 
(□), electron-AC-phonon (o) as well as electron-RI-phonon 
(T) scatterings as function of Fermi energy Ep. T = 200 K. 



In Fig. [T] we plot the total SRT and the contribu- 
tions from each individual scattering as function of tem- 
perature T. It is seen that the SRT changes little with 
temperature when T varies from 5 K to 100 K, and in- 
creases with increasing temperature when T > 100 K. 
The underlying physics is as follows. As said above, the 
inhomogeneous broadening does not change with tem- 
perature and electron density, therefore the temperature 
dependence of the SRT is determined by the momen- 
tum scattering: stronger momentum scattering leads to 
longer SRT in the strong scattering limit. ™ The electron- 
impurity scattering, which dominates the momentum 
scattering at low temperature, depends weakly on tem- 
perature. Thus the SRT varies very mildly with T. How- 
ever, the electron-phonon and electron-electron Coulomb 
scatterings both increase with temperature ^^'^Si'^i and 
become comparable to the impurity scattering at high 



temperature. This enhances the momentum scattering 
and hence increases the SRT. It is noted that the electron- 
electron Coulomb scattering, which is absent in the pre- 
vious investigations on spin relaxation in graphene^"— 
plays an important role in spin relaxation at high tem- 
perature. We also show that the electron-OP-phonon 
scattering is always negligible in the parameter regime of 
our investigation, which is consistent with the claim in 
the previous literaturei^'^ 

Then we turn to the electron-density dependence. In 
Fig. [21 the total SRT and the contributions from var- 
ious scatterings are plotted against the Fermi energy 
Ep at T ^ 200 K42. It is seen that the SRT decreases 
first rapidly and then slowly with increasing Ep. To 
understand the underlying physics, we first discuss the 
electron-density dependence of each individual scatter- 
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ing. The electron-impurity scattering decreases with 
due to the decrease of the cross section4i. The electron- 
electron Coulomb scattering also decreases with Ep in 
the degenerate regime due to the increase of the Pauli 
blocking»2£ The electron- AC-phonon scattering increases 
with increasing Ep since the matrix element (~ q) and 
the density of states (~ k) both increase with Ep ^ The 
electron-RI-phonon scattering varies slowly with Ep due 
to the competition of the decrease in the matrix element 
and the increase in the density of states4^ Under the joint 
effects of these factors, the behaviour in t — Ep curve is 
understood: the SRT first decreases rapidly with Ep due 
to the decrease of the electron-impurity scattering, and 
then decreases slowly since the increase of the electron- 
AC-phonon scattering partially compensates the decrease 
of the impurity scattering. 

The initial spin polarization dependence of the SRT is 
also investigated. In Fig. [31 we plot the SRT versus initial 
spin polarization P a.t T — 20 K and 200 K. It is seen 
that the SRT increases rapidly with the increase of the 
initial spin polarization. By comparing the calculation 
with and without the Coulomb HF term, one can see that 
the increase of the SRT originates from the Coulomb HF 
term. The underlying physics is similar to the previous 
studies in semiconductors>i^ the Coulomb HF term serves 
as an effective magnetic field along the z-axis, which is 
described by 

i?HF(k) = Yl ^kl'^kk'(4k'i - 4k'-l)/ (10) 
/ik' 

This effective magnetic field blocks the spin precession 
and slows down the spin relaxation. It is also shown 
that the SRT increases slower with P when temperature 
increases. This is because at high temperature the elec- 
trons are distributed to a wider range in k space, thus the 
effective magnetic field becomes smaller [see Eq. ([TU])] and 
the effect of the HF term is weakened?^ It is noted that in 
graphene there is a considerable increase of the SRT with 
the initial spin polarization, even at room temperature. 
In contrast, in semiconductors the electron system is in 
the nondegenerate regime at room temperature (as the 
Fermi energy in semiconductor is only tens of meV), and 
thus the effect of the HF term becomes insignificant 
This means that the HF effective field is more pronounced 
in graphene compared with semiconductors. Since the ef- 
fects of the HF term have been probed experimentally in 
semiconductors recentlyji^^— we expect they can be ob- 
served easily in graphene. 

We also study the effect of the in-plane electric field 
on spin relaxation. In Fig. 21 the SRT and the in-plane 
electric field £^|| are plotted against the hot-electron tem- 
perature Te for lattice temperature T = 300 K and ap- 
plied magnetic field B — 2 T. The hot-electron temper- 
ature is obtained by averaging the inverse of the slopes 
of g^ks = log(l//^k5 ~ 1) "^ith k varying along different 
directions. Here f^^^^ represents the hot-electron distri- 
bution function in the steady state, whose expression will 
be discussed in the next subsection. From Fig. [H it is 
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FIG. 3: (Color online) SRT as function of initial spin polar- 
ization P with (solid curves) and without (dashed curves) the 
Coulomb HF term at T = 20 K (•) and 300 K (□). 

seen that the SRT increases monotonically with the elec- 
tric field, which is very different from the complicated 
behavior in semiconductor quantum wellsJ^i^ In those 
systems, the high electric field induces two effects: (i) the 
drift of the electron distribution which enhances the inho- 
mogeneous broadening as more electrons are distributed 
at larger k and the SOC increases with k; (ii) the hot- 
electron effect which enhances the momentum scatter- 
ing. The former tends to decrease the SRT while the 
latter tends to increase. Therefore the electric field de- 
pendence of the SRT can be nonmonotonic J^i^i However, 
the spin-orbit field in graphene is independent on magni- 
tude of k. Thus Effect (i) on spin relaxation is marginal, 
and the electric field dependence of the SRT is mainly 
from the hot-electron effect. Consequently the SRT in- 
creases monotonically with £^|| . In addition, by compar- 
ing the calculation with and without the electron-electron 
Coulomb scattering at the same hot-electron tempera- 
ture, one finds that the contribution to the SRT from 
the electron-electron scattering increases with increasing 
hot-electron temperature, which is consistent with the 
lattice temperature dependence discussed above. 

B. Steady-state hot-electron distribution function 

Previous investigations have shown that in the Boltz- 
mann limit, the strong electron-electron Coulomb scat- 
tering can establish the steady-state hot-electron distri- 
bution function^ 

fX = {exp [(e^k - u • k - tis)/ikBTe)] + 1}"' , (11) 

where /is stands for the chemical potential of electrons 
with spin s and u is the drift velocity. Recently Bistritzer 
and MacDonald applied this distribution function to 
study the charge transport in graphene.'*'* However, 
whether the Coulomb scattering in graphene is strong 
enough to justify the validity of this distribution function 
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FIG. 4: (Color online) SRTs (•) and in-plane electric field iSy 
(□) as function of hot-electron temperature Te with (solid 
curves) and without (dashed curves) the electron-electron 
Coulomb scattering at lattice temperature T — 300 K and 
applied magnetic field in the Voigt configuration B = 2 T. 



remains unchecked. Since we can obtain the distribution 
function with the genuine Coulomb scattering explicitly 
computed, we check the validity of Eq. here. It 
is noted, if Eq. ^ is valid, g^ks = log(l//^L - 1) = 
(wF|k| — u-k— /is)//cBTe has a minimum around fc = and 
different slopes along different /c-directions. In Fig. [5ja), 
<?/jks is plotted against k varying along the direction of 
the electric field at E'y = 2 kV/cm and T = 300 K. One 
immediately finds that the hot-electron distribution func- 
tion from our calculation is very different from Eq. (1111) : 
the minimum of (7^ks is away from the point of fc = 
and the slopes are close to each other when k varying 
along opposite directions. Based on the above property, 
we propose the approximate expression of the computed 
hot-electron distribution function as: 



/uL = {exp [(e^k-u - f^s)/kBTe] + 1}' 



(12) 



Correspondingly, g^ks = (wplk - u| - ^s)/kBTe. To ex- 
amine our assumption, we also plot .g^^ks from Eq. (jl2p in 
Fig. [SJa) and find that the computed hot-electron distri- 
bution function is in reasonable agreement with Eq. (IT2t . 
In fact, for systems with parabolic energy dispersion, 
e.g., semiconductors in our previous investigations) ""^^i^" 
Eqs. (|lip and (|12p are equivalent. However, for graphene 
with linear dispersion, these two distribution functions 
are quite distinct. It is noted that Eq. (IT^ is just used 
to estimate the hot-electron temperature. The SRT and 
hot-electron distribution in this investigation are explic- 
itly computed from the KSBEs. 

In order to reveal the effect of Coulomb scattering to 
the steady-state hot-electron distribution, we introduce a 
dimensionless scaling coefficient x in front of the electron- 
electron Coulomb scattering, with x = 1 corresponding 
to the genuine case. We plot g^ks against k along the 
direction of the electric field with different scaling coeffi- 
cients X in Fig. [SJb). It is seen that with the increase of 
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FIG. 5: (Color online) (a) g^y^i = \og{l/f^^i - 1) from the 

KSBE computation (red dots) and from Eq. (|12p (blue solid 
curve) against k varying along the direction of the in-plane 
electric field at = 2 kV/cm and T = 300 K. (b) g^^i vs. k 

with the scaling coefficients x — 1 (red solid curve), 2 (green 
dotted curve) and 10 (blue dashed curve). 



X, i.e., the electron-electron Coulomb scattering, the dis- 
tribution function gets closer to Eq. pT|) . This indicates 
that in graphene the electron-electron Coulomb scatter- 
ing is not strong enough to establish the hot-electron dis- 
tribution Eq. pT|). 



IV. CONCLUSION AND DISCUSSION 

In conclusion, we have investigated the spin relax- 
ation in graphene from the microscopic KBSE approach, 
where all the relevant scatterings, especially the electron- 
electron Coulomb scattering, are explicitly included. We 
show that the SRT remains almost unchanged with in- 
creasing T at low temperature because the electron- 
impurity scattering, which dominates the momentum 
scattering, varies little with temperature. Nevertheless, 
the SRT increases with T at high temperature because 
the electron-electron and electron-phonon scatterings be- 
come comparable to the electron-impurity scattering and 
both scatterings increase with increasing T. We also 
show that the electron-electron Coulomb scattering plays 
an important role in spin relaxation at high temperature. 
It is also seen that the SRT first decreases rapidly with 
the increase of Ep due to the decrease of the electron- 
impurity scattering, and then decreases mildly with Ep 
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since the increase of the electron-AC-phonon scatter- 
ing partially counteracts the decrease of the electron- 
impurity scattering. We also predict a pronounced in- 
crease of the SRT at high spin polarization in the whole 
temperature regime of our investigation. The underlying 
physics is that the Coulomb HF term serves as an ef- 
fective longitudinal magnetic field which blocks the spin 
precession and suppresses the spin relaxation. The ef- 
fect of the in-plane electric field on spin relaxation is also 
investigated. It is shown that the SRT increases with 
the in-plane electric field due to the hot-electron effect. 
Moreover, we show that the electron-electron Coulomb 
scattering in graphene is not strong enough to establish 
the usual steady-state hot-electron distribution in the lit- 
erature and an approximate one is suggested based on our 
computation. 

Now we address the effect of ripples on the spin relax- 
ation. For graphene samples with an undulating sur- 
face, i.e. ripples," an additional Rashba-type SOC, 
whose expression is the same as the electric field induced 
SOC^ appears. The SOC coefficient due to the curva- 
ture effect is estimated to be acurv~0.02 meV (0.2 K))^ 
which is about two orders of magnitude larger than 
ttR = 1.5 X 10"'^ meV used in our calculations. For the 
well-known relation l/r oc for the DP mechanism;^ 
one expects that the SRT is shortened by four orders of 
magnitude. Nevertheless, the additional SOC docs not 
change the dependences of the SRT on the temperature 
and sample parameters, as well as the importance of the 
electron-electron Coulomb scattering. 

Finally, it is noted that even after considering the 
SOC enhanced by the curvature effect, the SRTs 
from our calculation are still two orders of magnitude 
larger than those in recent spin transport experimental 
measurements <^">22, We stress that the range of impurity 
(adatom) density of the graphene samples we discuss is 



different from that in recent experiments, and therefore 
the dominant spin relaxation mechanism is also quite dif- 
ferent. As mentioned above, the mechanism most likely 
limiting the SRTs in the recent experiments is the lo- 
cal spin-orbit field induced by the randomly distributed 
adatoms. The experimental and theoretical works22ii& 
showed that the SOC strength from the adatoms can 
reach 10 meV, which is about three orders of magni- 
tude larger than the strength used in our calculations. 
The study on the effect of the adatoms on spin relax- 
ation is beyond the scope of this investigation. It is 
further noted that the spin-related experiment in clean 
graphene is still missing, we expect that the effects pre- 
sented in this manuscript can be confirmed by the fu- 
ture experiments in relatively cleaner graphene samples. 
A possible method to obtain the graphene sample with 
higher mobility is the epitaxial growth of graphene on SiC 
substrate42r— Nevertheless in this system, some other 
factors must be taken into account, e.g., the scattering 
arising from the interfacial states4^ The study in this 
system can be the future extension of this investigation. 
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Appendix: Expression of the interaction Hamiltonian 

The electron-impurity and electron-electron Coulomb 
interaction Hamiltonian can be written as 



^ei — / , ^k.k-g^l 



kk-q^ ^ks ^k-qs' 



(A.l) 



ks ^\<i's' '^k'+qs'^^k-qs 



yuu' rp^uy' rp^' u' y fi'i^'' y' v yu' 



, ^ TT-llUu' rpflViy' rpfl'v' 

9 2^ ^k,k-q^kk-q^k'k' 



+q'^ks *^k's' '^k'+gs'^k-qs- 



(A.2) 



In these equations, Rj stands for the position of jth im- 

purity; T^-^q = ry'r^U ^Ck-q ^^^k':k-qe-''^ is 
the electron-impurity interaction matrix element. Here 
= 1 is the charge number of the impurity; the ef- 
fective distance d of the impurity layer to the graphene 
sheet is chosen to be 0.4 nm,ii>^>^-^ ^^k'^k-q denotes the 



screened Coulomb potential where the screening is calcu- 
lated under the random phase approximationj^LMi^^i^ 



k,k-q 



(0) 



yvV. ^/iz.''k— q 



(A.3) 
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where Vc 



(0) 



2'KVYrs/q is the two-dimensional bare 



Coulomb potential with rg = 0.8i^"— As pointed out in 
Refs. [13 andim, such small ensures the validity of the 
random phase approximation. n(q, to) is given \yy^^— 



n(q,^)= E 



^rp^VU |2 

Kkk-ql 



/xi-'i/'ks 



^/j.^'k ^/j.t^'k— q 



(A.4) 



It is noted that the interband contribution in screening 
cannot be neglected even in the n-doped case. 

The electron-phonon interaction Hamiltonian takes the 
form (only the terms relevant to the electron band pre- 
sented) 



Hp 



E 

fifi' ,u=i 

kqs 



(A.5) 



Here q (lA.q^) is the annihilation (creation) op- 
erator and M^y^ M'k-q stands for the matrix element 
of the electron-phonon interaction with A being the 
phonon branch index. For the electron-AC-phonon 
scattering, the phonon energy spectrum uj^'~^ — Vp^q 



and \M^^y^^ 



rv-^kk-q'^A'M'' where Vph 



2 X 10^ cm/s is acoustic phonon velocity and Pm = 
7.6 X 10^* g/cm^ denotes graphene mass densityi^i^ 
The value of deformation potential D is still in de- 
bate. The range of D obtained via various theoretical 
and experimental methods is from 4.5 to 30 eVj ^^'^"" — 
Here we choose a moderate value D — 19 eV.— >^i^ 
For the electron-RI-phonon scattering, \Mj^^^^,^_^\'^ = 

kk-q^/i/i' where Qs — ^rghp is the Thomas- 



q— 



/k 



Fermi screening length, 



57 



1.42 A is the C-C bond 



distance. For Si02 substrate, the energy of the re- 
mote phonon modes are ujf^ — 59 meV and lu^^ = 
155 meV; the corresponding dimensionless coupling pa- 
rameters are gi — 5.4 x 10~^ and 32 — 3.5 x 10~^, 
respectively!^ The matrix element of the electron-OP- 
phonon scattering is described by l-^/^f^'k-qP = 2^1^- 
For the longitudinal and transverse optical phonons near 
r point, which cause the intravelley scattering, AY~' = 



(Dl)[l - cos(0k + ^k-q - 20q)]V, ^™ = (^r>[l + 
cos(6ik -I- 6'k-q - 26lq)](5^^/ with (Dl) = 45.60 eV^/k^ 
and ujr = 196.0 meV; whereas for the transverse optical 
phonon near K{K') point which causes the intervalley 
scattering, Aj^ = - cos(6'k - ffk-q)]^^,-;^' with 



92.05 eVVA^ and ujk = 161.2 meV.-^i^ 
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